Angular velocity distribution of a granular planar rotator in a thermalized bath 
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The kinetics of a granular planar rotator with a fixed center undergoing inelastic collisions with 
bath particles is analyzed both numerically and analytically by means of the Boltzmann equation. 
The angular velocity distribution evolves from quasi-gaussian in the Brownian limit to an algebraic 
decay in the limit of an infinitely light particle. In addition, we compare this model with a planar 
rotator with a free center. We propose experimental tests that might confirm the predicted behaviors. 
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I. INTRODUCTION 

When macroscopic particles undergo inelastic colli- 
sions the total kinetic energy decreases with time. If 
an external source of energy, such as a vibrating bot- 
tom wall, is present, the system may reach a station- 
ary state. Despite similarities with equilibrium sys- 
tems, however, equilibrium statistical mechanical con- 
cepts cannot be applied 0, 0, El For instance, there is 
no equipartition between different species in polydisperse 
systems 0, H, IE S 0: an d velocity distributions are, in 
general, non- gaussian [ii fToL Hll| . 

In dilute systems, most collisions only involve two par- 
ticles, and consequently, a theoretical description of the 
dynamics has been proposed starting from the Boltz- 
mann or Enskog eauation[3.IT^|. 

The inelastic hard sphere model is a paradigm for gran- 
ular gases in the same way that the elastic hard sphere 
model is for equilibrium fluidsQ- While both account 
for the exclusion effect, energy is lost at each collision in 
the former. For dissipative systems, recent studies have 
shown the relevance of the granular temperature, even if 
the absence of equipartition yields a number of temper- 
atures increasing with the polydispersity or the number 
of degrees of freedom of each particle. 

Few studies have addressed the effect of particle shape 
on the properties of granular gases. In three dimen- 
sions, Aspelmeier et al0 showed that the rotational and 
translational temperatures are different in the free cool- 
ing state of inelastic hard needles. Anisotropic trac- 
ers (needle and spherocylinder) in a thermalized bath 
also displa y n on-equipartition between different degrees 
of freedom [13. Il4j . The purpose of this paper is to build 
a simple model which captures the specific features due 
to the shape of the particle and the presence of an irre- 
versible microscopic dynamics. 

We present here an investigation of the stationary ki- 
netics of a planar rotator with a fixed center that un- 
dergoes inelastic collisions with the bath particles. With 
these assumptions, the kinetics is described by a linear 
Boltzmann equation. We show that, when the tracer is 
much heavier than the bath particle, the angular veloc- 



ity distribution function is quasi-gaussian (the Brownian 
limit) whereas it exhibits an algebraic decay in the oppo- 
site limit of an infinitely light granular particleflij. For 
all intermediate cases, there is no simple scaling regime 
and deviations from gaussian behavior arc captured by 
analyzing the fourth moment of the distribution function. 

The paper is organized as follows: the model, its me- 
chanical properties and the Boltzmann equation are pre- 
sented in section [H] The asymptotic solution of the 
Boltzmann equation in the Brownian limit is presented 
in section 11111 In section I1VI analytical results are de- 
rived for a zero-mass limit and for different coefficients 
of restitution and intermediate cases are considered in 
section In section IVII we compare the fixed rotator 
with one whose center is free, both within the gaussian 
approximation. Finally, the conclusion discusses possible 
experimental tests of our theoretical results. 

II. THE PLANAR ROTATOR 
A. Definition and mechanical properties 

The model consists of a two-dimensional, infinitely 
thin needle of mass M, length L and moment of iner- 
tia / = ML 2 /12 immersed in a bath of point particles, 
each of mass m. The needle has a fixed center of mass, 
but can rotate freely around its center. It undergoes in- 
stantaneous and inelastic collisions with the surrounding 
bath particles. The motion of the planar rotator can be 
described with the angle between a unit vector u collincar 
to the axis of the needle and the x-axis. 

The rate of change of the orientation u = ujuj_ is equal 
to the angular velocity uj £] — oo, +oo[ times a unit vector 
ui perpendicular to u. 

The angular velocity of the rotator changes at each 
binary collision with a bath particle. The position of 
the point of impact along the needle axis is denoted by 
Au. Obviously, a condition for collision is |A| < L/2 (see 
FigfjJ. The relative velocity V at the point of impact 
equals 

V = v Aii = v — Xuju± (1) 



2 



where v denotes the velocity of the bath particle. 

Due to the dissipative nature of the collision, the rela- 
tive velocity changes according to the collision law 



-aV± 



V\ 



(2) 
(3) 



where < a < 1 is the normal restitution coefficient and 
where the indices _!_ and 1 1 indicate the perpendicular and 
parallel components of any vector relative to the needle 
axis, respectively. When a — f, one recovers an elastic 
collision rule. For the sake of simplicity, the tangential 
component of the velocity is unchanged during the colli- 
sion. 

Since each collision conserves the total angular momen- 
tum we have that 



Iuj* + \mv*i = Iu) + Xmvj 



(4) 



where post-collisional quantities are denote with a star. 

By combining Eqs. Q-Q, the post-collisional bath 
particle velocity is given by 



v ± = vj_ 



1(1 + a)V± 
I + mX 2 



(5) 



whereas the corresponding post-collisional angular veloc- 
ity is 



* i 

LO = LU + 



(f + a)V±mX 
I + mX 2 ■ 



(6) 



The inverse transformation (giving the pre-collisional 
quantities denoted by a double star) is obtained by sub- 
stituting a by a -1 and the starred quantities by double- 
starred quantities. 



B. Homogeneous Boltzmann equation 

At low density, one assumes that the needle influences 
weakly the local density of the bath and, consequently, 
that the system remains homogeneous. After a transient 
time (not considered here) the kinetics of the needle be- 
comes stationary and can be described by the stationary 
Boltzmann equation. This expresses the invariance of the 
rotator angular velocity distribution function, F(uj), re- 
sulting from a balance between collisional gain and loss 
terms: 



j-L/2 j 

/ dX dv\v± — Xlu\ 

J-L/2 J 



F(w**)$ B (v**) 



= 0,s 



-F(w)* B (v) 
(7) 



where the pre-collisionnal velocities v** and co** are given 
by the right hand sides of equations (JjJJ and ffifl. respec- 
tively, with a replaced by oT 1 . $b(v) is the time inde- 
pendent bath velocity distribution. 




FIG. 1: Illustration of a collision between the planar rotator 
and a bath particle: u and u x are unit vectors parallel and 
perpendicular to the axis of the rotator. The element of the 
needle at Au moves with linear velocity Xlou . 



The integration over the parallel velocity component 
Wii in Eq.Q can be readily carried out since F does not 
depend of this variable. By changing of the integration 
variable we can rewrite Eq.Q as 



L/2 



-L/2 



dX / dv±\v±\ 



F ( w + v± 



(I + a)mX 



Xlu + v± 



(amX 2 — I) 



f mX 2 

with 4>b(v) = J du||$s(|v| 



I + mX 2 
F(lu)4>b(v± + Xui) 



0. 

(8) 



III. BROWNIAN LIMIT 

An exact solution of Eq. (JSJ cannot be obtained in gen- 
eral. When the mass of the planar rotator is much larger 
than the mass of the bath particle, however, one expects 
that the deviation from Maxwellian behavior is weak. It 
turns out (see Appendix ^) that exploring the regime 
corresponding to Brownian motion is equivalent to the 
analysis of the small A expansion of the collision term 
(8). We thus perform perform a pcrturbative expansion 
of the integrand of Eq.JSJ in terms of A. We denote 



G(v_l,uj,X) 



F 



(j)B [ XtU + V± 



(amX 2 — I) 
I + mX 2 



(I + a)mX 
+ V± I + mX 2 

- F(lu)4>b{v±_ + Xtu) 



(9) 
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with the property that G(v±,uj, 0) = 0. To go fur- 
ther, we assume in the rest of this section that the bath 
distribution, </>#(«), is a Maxwcllian: 



4>b{v) = 4>m{v) = 



2ttT 



exp(-mv 2 /2T) 



(10) 



The first derivative of G(v±,uj, A) with respect to A at 
A = gives the differential equation 



^dF(uj) „, ^ujI 
(l + a)-^+2F(u) — 



whose solution is 

F(u>) oc exp 



Iuj 1 



(l + a)T 



(11) 



(12) 



By taking the second derivative, one obtains the dif- 
ferential equation 



(1 + a) 



d 2 F{uj) 
duj 2 



21 dF(oj) 21 . , 

— uj — H F(uj) 

T du T 







(13) 



whose solution is also given by Ea. H12|) . 

The third-order expansion gives a polynomial in v± 



m 3 v 3 L 

where 
H 3 (F(u>)) 



H 3 (F(uj)) 



3m 2 vi 



■Hi (F(w)) = 



(14) 



(1 + a) 



I 



3 d 3 F(w) 
duj 3 
Iuj 2 



3(l + a) : 



Jw d 2 F(uj) 



T duj 2 



3- (1 + a) 

TV T 



2I 2 uj 



dF(w) 
du 



Iuj< 



.3(1 + a) + — ] F(w) 



(15) 



and 



+ 



H 1 (F{oj)) = (l + a)(2 + —)-±J. 

(1 + a) + 2^-) F(u) 



Iuj 
~~T 



T 



(16) 



The solution of H 3 = is given by Eq. l|12|l again, but 
the solution of H\ = is 



F(oj) 



Iuj 2 
T 



1 — o 
1+a 



exp 



-Iuj 2 



(l + a)T 



(17) 



If Ea. lfTTjl if different from Ea. (jT2*)l . it is interesting to 
note that this solution is the same Maxwellian times a 
slow decreasing function (If a = 1, Ea. ll7|) is identical to 
Ea. (|12f> ). The fourth derivative of G(v±,uj, A) has been 
calculated with the software Maple, and it occurs that 
the differential equation associated with the term has 
a Maxwcllian solution but the solution of the differential 
equations associated with the v±-tevm is given by the 



(i.tii 
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FIG. 2: (a) Log-linear plot F(uj) for M/m = 10 and 
with different values of the coefficient of restitution a = 
0,0.2,0.4,0.6,0.8,0.9 (from top to bottom, center). (b) Log- 
linear plot F(uj)exp(luj 2 /((1 + a)T)) for M/m = 10. The 
insets display the same curves on a larger range of angular 
velocities. 



Maxwellian solution multiplied by a slow varying func- 
tion. Therefore, we conjecture that the complete solution 
is given by Ea. H12(l times a sub-dominant term. 

In order to check this assumption, we have solved nu- 
merically the Boltzmann equation. Since Eq. JSJ) is linear 
and the distribution F(uj) is a one-variable function, we 
have used an iterative method that is very efficient and 
provides a much more accurate solution^] than a DSMC 
method[l3|. 

By using Eq.ljSJl, the numerical resolution is obtained 
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by iterating the following equation 

^" +1 >H =C(u;) [dxf dv\v\F™ ( w + v [l + a) ™ X 
J J V I + mX 2 

(am\ 2 — I) 



<f>M Aw + u- 



wherc 



C(w) 



/ + m\ 2 



dX / — Xuj\(/)m{v) 



(18) 



(19) 



C(w) is an explicit function when the bath particle distri- 
bution <Pm{v) is a Maxwcllian. The velocity distribution 
is sampled on a one-dimensional grid with 1000 points. 
The integrations over A and w are performed with Simp- 
son's rule with 100 and 1500 points, respectively. For 
values of the velocities which do not match to the lattice 
points a linear interpolation is performed. The initial 
distribution is taken as the Maxwellian, Ea. (|12|) . Except 
when the mass of the planar rotator is extremely small, 
the method converges very rapidly. 

Fig[2Ji. displays the angular velocity distribution F(ui) 
for a mass ratio M/m = 10 for different values of the 
coefficient of restitution a. The deviations from the 
Maxwellian, shown in Figl^b, increase with decreasing 
coefficient of restitution, but compared to the function 
F(uj), they vary weakly with the angular velocity. The 
insets show that the iterative method allows the tails of 
the distribution function to be obtained with high preci- 
sion (unlike with the DSMC method) . Fig|3] shows the 
angular velocity distributions for a mass ratio equal to 
one. The deviations from the gaussian behavior are more 
pronounced than those for a mass ratio of 10, but they 
are still negligible compared to the the leading gaussian 
term. 



IV. ZERO-MASS LIMIT 

When both the mass of the needle and the coefficient 
of restitution are equal to zero, the Boltzmann equation, 
Eq.lJSJl, can be solved exactly. Let us review the special 
characteristics of this limiting case. First, because of the 
zero mass of the tracer particle, the velocity of the bath 
particles never changes as the result of a collisions. In ad- 
dition, the angular velocity of the needle has no memory 
of its pre-collisional value. This quantity is reset after 
each collision, acquiring instantaneously the value v±/X. 
By introducing the variable v± — (Xto)y, Eq.JSJ) can be 
written as 



L/2 



d\\' 



L/2 

B [ Xui I 1 + y 



J dy\y\F(w U + V^j 
(amX 2 - I) 



(1 + a)mA 5 



(A 2 



I + mX 2 



L/2 



= F(u) I dXX 2 I dy\y - l\</> B (Xu}y) 

-L/2 



(20) 
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FIG. 3: Same caption as Fig.2 for for M/m = 1 



When 7 = and a — 0, this simplifies to 



dXX 2 / dy\y\F(w(l + y))(f> B (Xu) 



L/2 



-L/2 

,L/2 , 

= F(w) / cfAA 2 / dy\y + l\4> B {Xwy) (21) 

J-L/2 J 

The exact solution F(u>) of Ea. H21(l is obtained by in- 
tegrating the bath distribution (f> B (Xu>) weighted by the 
position A of the point of impact: 







[ L/2 dX 




J-L/2 





(22) 



5 



o.oi 



3 



0.0001 - 



le-06 - 




FIG. 4: Log-linear plot F(lj) of the "zero mass" needle with 
uj c — 2. The upper inset displays the log-log plot showing the 
power-law decay for sufficiently large u> values (the vertical 
dashed line corresponds to uj = ui c ). The lower inset is a 
log-log plot showing the gaussian behavior for uj < u c (the 
dashed curve corresponds to the gaussian approximation) 



Note that this solution is independent of specific as- 
sumptions about the distribution of the bath particles. 
If we make the weak assumption that the second mo- 
ment of the distribution 4>b(v) is finite, i.e. that the 
bath is characterized by a finite granular temperature, 
it can be shown that F(ui) decays algebraically as oj~ 3 . 
The long tail of F(uj) arises from collisions near the cen- 
ter of the rotator that result in large angular velocities. 
It is worth noting that solutions of the Bolzmann equa- 
tion with power-law decay exist for isotro pic p articles, 
but not with a thermalized bath of particles [LSI Il9l| . 

While the granular temperature of the bath particles is 
finite, that of the tracer is not well defined since its zero 
mass implies an infinite mean squared angular velocity. 
This difficulty is removed when the granular particle has 
a small, but finite mass (see below) 

An explicit expression of the angular velocity distribu- 
tion function F(uj) can be obtained when the bath par- 
ticles have a Maxwellian velocity distribution, 4>b{v) = 
^Jm/2nT exp(— mv 2 /2T) . In this case 



only decreased by an order of magnitude compared to 
the maximum, F(0). 

An obvious question is whether the behavior just de- 
scribed persists when either a or M/m is different from 
zero. Is the power law regime sustained in these cases? 

We first consider the case where the mass ratio is main- 
tained at zero, but the coefficient of restitution is allowed 
to take all values between and 1. For a = 1 (clastic 
collisions) the solution of the Boltzmann equation is the 
expected Maxwell distribution: 



2irT 



exp 



2T 



(24) 



but the limit / — > does not yield a probability distri- 
bution. 

For < a < 1 and / = 0, we were unable to obtain 
an analytical solution and we investigated the behavior 
numerically with the method described above. Figure 
[3] shows the logarithm of the distribution function ver- 
sus the angular velocity. As in the case a = 0, two 
distinct regimes characterize F(tu): a scaling regime for 
uj > uj c (pt), where uj c (a) is a cutoff which increases with 
a and a gaussian behavior at low frequencies. 

FigEl shows that, in the power-law regime, while the 
amplitude of F(u>) decreases when a increases, the ex- 
ponent of the power-law is independent of a. It is pos- 
sible to obtain this result analytically by means of an 
asymptotic analysis of the Boltzmann equation, Ea. (|21|l . 
Details of the calculation are given in Appendix [BJ The 
final result, for M = and < a < 1, is 



4T(1 + a) 
mL 2 (l — a)uj 3 



(25) 



Unlike the case where a = 0, the rotator has a memory 
of its previous angular velocity after a collision with a 
bath particle, uj* = auu + (1 + a)v±/X. This memory 
effect can be very small when a bath particle collides 
near the center of the rotator (small A), which explains 
why, for large angular velocities, the distribution function 
behaves similarly to the case a = 0. 

The inset of FigEl displays H(uj) = F(u) ^^^ w 3 
versus uj, showing that the asymptotic behavior is rapidly 
reached and that the numerical results agree accurately 
with Eq. Qgfy . 



INTERMEDIATE CASES 



F(uj) = - 




where uj c = ^8T/(mL 2 ) is a crossover frequency be- 
tween two regimes. For uj > uj Ci F{u>) has a power-law 
behavior and for uj < oj C} F(uj) ~ cxp(— 3ui 2 /(5uj 2 )) is 
gaussian (see Fig0J. It is important to note that when 
the power-law regime begins, the amplitude of F(uj) has 



We have considered in the preceding sections the two 
limiting cases of a large mass of the planar rotator (Brow- 
nian limit), where the angular velocity distribution func- 
tion displays a gaussian-like behavior and a zero mass 
case where, surprisingly, a solution of the Boltzmann 
equation exists with a power-law decay of the distribution 
function. In this section, we investigate the intermediate 
case that corresponds to most physical situations. 
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FIG. 5: Log-Log plot of F(uj) of the "zero mass" needle 
with different values of the coefficient of restitution a — 
0, 0.2, 0.4, 0.4, 0.8 from top to bottom. The inset shows the in- 
dependence of the rescaled function H{u) = F(u) "^{v+a)^ 
for ujc — 2. 



FIG. 6: Log-Log plot of F(uu) for a needle with a mass M/m = 
0.005,0.01,0.1 (dashed curve corresponds to M = 0.0). The 
inset shows the crossover between the power law and the 
gaussian-like asymptotic 



We consider a planar rotator with a small but finite 
mass < M <C m. In Ea. i|20[l . since the integrand is 
an even function of A, integration can be restricted to 
the positive abscissa. Moreover, the integration can be 
divided into two parts: 



L/2 



dX 



dX 



L/2 



dX 



(26) 



where e > > y (but e being always small compared 
with L/2). For small angular velocities co, the contri- 
bution of the first integral vanishes, and performing a 
first-order expansion of the arguments of the integrand, 
one obtains 



L/2 



d\F{u(l + y(l + a)))<l>B(\(l + ay)) (27) 



which leads to the Boltzmann equation of the massless 
particle when e — > 0. 

The existence of a finite mass allows for restoring 
a finite granular temperature. For small mass M, as 
shown above, the angular distribution function F(u>) has 

a power decay regime truncated at an upper ui c where 
the gaussian begins. The granular temperature given by 
the product of the momentum of inertia times the second 
moment of the distribution function is dominated by the 



integral in the power-law regime. 



T = J dujF(uj)Iuj 2 
7^~71n(I) 

,(2) U) 



(28) 



which means that the granular temperature of small- 
mass planar rotator goes to zero when the mass decreases 
(even if the quadratic average of the angular velocity di- 
verges logarithmically) . 

FigEI shows the distribution function for three small 
mass ratios M/m = 0.005,0.01,0.1 with a = 0. As 
expected the low frequency distribution is well approxi- 
mated by the massless distribution function. The inset 
shows that the range of the power-law decay decreases 
as the mass of the planar rotator increases. For large 
angular velocities, the distribution function resumes the 
gaussian behavior, F(w) ~ exp(— Icu 2 /((1 + a)T)), irre- 
spective of the needle-to-bath particle mass ratio. 

In summary, the velocity distribution keeps memory of 
the massless solution up to second crossover angular ve- 
locity lu ( c 2) ~ ^/2T/1. If the mass of the planar rotator is 
sufficiently small, one observes three successive regimes: 
first, a gaussian decay, second a power-law and finally a 
gaussian-like decay (sub-dominant terms arc present). 

When the masses of the planar rotator and a bath par- 
ticle are comparable, the two crossover angular velocities 
merge and the power-law regime disappears. However, 
F{lo) may still deviate significantly from a gaussian. In 
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FIG. 7: Deviations of the angular velocity distribution func- 
tion from a gaussian for M/m = 0.1, 1, 2, 10, top to bottom. 



order to show this we introduce the quantity 



V 



i(u 2 f 



(29) 



which is zero for a gaussian distribution. FigQdisplays ?/ 
as a function of the coefficient of restitution a for different 
masses M/m = 0.1, 1, 2, 10. 

For the spherical tracer particle, the distribution func- 
tion is a pure gaussian when the bath of particles is a 
gaussianjfj. Deviations from the gaussian occur, how- 
ever, in a mixture composed of granular particles of 
finite density immersed in a thermostat providing the 
energy via elastic collisions when this energy is redis- 
tributed within the granular component through inelastic 
encounters 1 16j|. For anisotropic particles, our calculations 
show that deviations are present even in the limiting case 
of infinite dilution. 



VI. MODELS WITH FREE OR FIXED CENTER 

The model of granular planar rotator with a free center 
in a thermalized bath has been previously investigated by 
two of us|l3(. In this model, the tracer particle has, in 
addition to the rotational one, two translational degrees 
of freedom. 

By using an approximate theory, as well as numerical 
simulations of the Boltzmann equation, it was shown that 
the translational and rotational granular temperatures 
are both smaller than the bath temperature and also dif- 
ferent from each other. In addition, the translational and 
rotational degrees of freedom are correlated [20j. 



In this section, we compare the rotational granular 
temperatures for the two models, all parameters being 
the same (mass, coefficient of restitution, bath tempera- 
ture,...). 

In order to obtain an analytical expression for the ro- 
tational temperature, we use a method originally pro- 
posed by Zippelius and colleagues that consists of calcu- 
lating the second moment of the angular velocity of the 
Boltzmann eauation[7l Il3l l2ll | . In a stationary state, this 
quantity is constant, and by using an gaussian ansatz for 
F(u) 



-lu? 

F{uj) oc exp ( - 1 = r 



(30) 



one obtains a closed equation for the granular tempera- 
ture T as a function of microscopic quantities: 



A T 2 

dx—x 
T 



-kx 2 



1 + kx' 



dx 



x 2 (l + ^kx 2 fl 2 



(1 + kx 1 



where 



k 



mL 2 
41 



(31) 



(32) 



Details of the calculation are given in Appendix [U] 
Eq. 1)3 1J) is an implicit equation for T, but, for a given 
value of a and of the mass ratio, it can be solved with 
standard numerical methods. 

Fig [5] shows the variation of the rotational tempera- 
ture with the normal coefficient of restitution for a mass 
ratio of 1 for a fixed (solid curve) and free (dashed curve) 
planar rotator. The temperature is always higher when 
the center is fixed except in the case of clastic collisions, 
a = 1. The circles correspond to the "exact" tempera- 
tures obtained by computing the second moment of the 
distribution function F(lo), which shows that the above 
method provides accurate approximate results for esti- 
mating the granular temperatures. 

We also consider the variation of the granular temper- 
ature with the mass ratio for a given value of the resti- 
tution coefficient: See Fig EI It is easy to verify that in 
the Brownian limit, i.e. when the mass of the rotator is 
much larger of the particle mass, the granular tempera- 
ture goes to a value, (1 + a)/2, that is independent of 
the mass ratio and, probably, the shape of particle. As 
the mass of the needle decreases, the difference between 
the rotational temperatures of the free and fixed planar 
rotator increases. For elastic systems, the rotational tem- 
peratures remain identical for the two different situations 
(free or fixed center of mass), but significant differences 
occur for a granular particle. 

This phenomenon is more pronounced with varying 
mass ratio than with varying coefficient of restitution. 
Therefore, by monitoring the rotational motion of a gran- 
ular needle in a bath of significantly heavier particles 
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FIG. 8: Effective granular temperature T as a function of 
the restitution of coefficient a for a fixed (full curve) and 
free rotator (dashed curve) in a bath of point particles. Cir- 
cles correspond to the exact temperatures by computing the 
second moment of the distribution function obtained by the 
numerical resolution. The mass ratio is M/m = 1 

in two successive experiments (free and fixed center of 
mass), it should be possible to observe the absence of 
cquipartition for a single particle. 

VII. CONCLUSION 

We have shown that the stationary angular velocity 
distribution of a planar rotator with a fixed center that 
collides inclastically with particles in a thcrmalizcd bath 
displays a variety of behavior as the mass of the rota- 
tor is changed. Starting from a quasi-gaussian regime 
when the rotator is much heavier than the bath parti- 
cle, F(u>) shows significant deviations from the gaussian 
when the mass of the rotator is comparable to that of a 
bath particle. As the rotator mass decreases further, an 
intermediate power-regime also appears. 

We believe that these features are not specific to this 
simple model, but are generic for all granular systems 
containing non-spherical particles. Furthermore, we ex- 
pect that the non-gaussian character that is present in 
the idealized configuration of a thcrmalizcd bath will be 
amplified in in the presence of a granular (non-thermal) 
bath. Several experiments on intensely vibrated gran- 
ular systems using high speed photography [H, 122^ . im- 
age analysis and particle trackingjj, |23j have shown 
that many quantities, including granular temperature 
and velocity profiles, can be precisely measured. As re- 
cent experimental studies have demonstrated, some of 
the same techniques can also be applied to granular 



FIG. 9: Effective granular temperature T as a function of 
the mass ratio M/m for a fixed (full curve) and free rotator 
(dashed curve) in a bath of point particles. The coefficient of 
restitution is a — 0.9. 

We believe that the significant 
difference between the granular temperatures of a free 
and fixed rotator that we have identified (Fig|5J should 
be detectable using available experimental methods. 
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APPENDIX A: DIMENSIONLESS BOLTZMANN 
EQUATION 

In order to consider the Brownian motion regime, it 
is convenient to introduce dimcnsionless variables f2 = 
\J I /ksTuj and u± = y m/ ksTv±. 

In terms of these variables the integrand in the Boltz- 
mann equation (8) takes the form 

-F(n)<f> B (u± + efi) (Al) 

where e = y/mX 2 /! = ^jYlm/MX/L. 

It is thus clear by inspection that the expansion of 
the dimcnsionless collision term (IA1|) in powers of e is 
equivalent to the expansion of the collision term (8) in 
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powers of A. Notice that keeping the variable f2 fixed 
when exploring the region of m <C M corresponds to 
the Brownian motion asymptotics, as then the rotational 
energy Iuj 2 is maintained at a fixed ratio with the thermal 
energy k B T. 



APPENDIX B: ASYMPTOTIC BEHAVIOR OF 
THE ANGULAR VELOCITY DISTRIBUTION IN 
THE ZERO MASS LIMIT 

Performing the change of variable u = oj(y + 1) in the 
left-hand side of Ea. i|2U[l and u = toy in the right-hand 
side, one obtains 



L/2 



L/2 



dXX 2 I du\u — uj\F ((1 + a)u — aoj) 
4>b (A((l - a)uj + au)) 

,L/2 , 

= F(tu) / dXX 2 / du\u-Lu\(f> B (Xu) 

J-L/2 J 

(Bl) 



When \ui\ — > oo, one has 



du\u — uj\(f>B(Xu) 



X 



(B2) 



By taking the Fourier transform of Eq. i|B4(l and by using 
Eq. (IB6() , the asymptotic form of the Boltzmann equation 
becomes 

Hk) = F(ak) J L ^dX\X\^ B (^«) (B7) 



The integral of the right-hand side of Eq. (|B7|) can explic- 
itly performed 



L/2 



L/2 



dX\X\(t)l 



1 + a 



( (l + a)¥T 
a/iexp -2 — 

(l + a) 2 k 2 T s 
mL 2 



exp —2 



2 (l + a)yT E ( 12 (l + *Vk 2 T 



mL- 



mL 2 
(B8) 



where Ei(l, x) is the exponential integral. For small val- 
ues of k, Ea. (|B8|) behaves as 



L/2 



.1 + a 



(l + a) 2 T 



dWfoi^-^q) ^ 1 + 4 V ^ ' fc 2 ln(|fc|) (B9) 



-L/2 



Inserting Eq.||B9J) in Ea. (|B7p yields 



and therefore, the right-hand side of Eq.jBTJ) becomes 

2 



(B3) 



Performing a similar analysis for the right-hand side of 
Ea. (IBl|) . one gets the asymptotic relation of the Boltz- 
mann equation (for the zero-mass limit) 



L 



L/2 



L/2 



dXX J duF((l + a)u - atu) 

4>b (A((l - a)u + au)) 

(B4) 



Let us introduce the Fourier transforms of the distri- 
bution functions F(cu) and 4>b{u) (with the convention 
F(k) = J dujF(uj)e- tkr ). The right-hand-side of Ea. (jB4|) 
can be expressed as 



dk f dq 



2tt / 2tt 



L/2 



L/2 



dXX 2 F{k)4>B{q) 



d ue ik l( 1 + a ) u - au1 } e i i x [( 1 - a ') ul + au ] (B5) 



By using the property J dxe lax = 2nS(a), integration 
over v in Ea. (|B5p can be carried out and one obtains 



/ * f L/2 dXX 2 P ( ^\ kM e H^) (B6) 
J 2k J-l/2 \l + aj 1 + a V 7 



F(k) 



F(afc)(l + 4 



{l + a) 2 T u2 



mL 2 



k 2 \n(\k\) (BIO) 



Iterating Ea. l|B10|) and by using that F(k) = 1, one ob- 
tains that 

F{k) = l + \±^^k 2 H\k\) (Bll) 
1 — a mL^ 

The inverse Fourier transform of Ea. (|Bll|l leads to 
Eq.® 



APPENDIX C: GRANULAR ROTATIONAL 
TEMPERATURE OF THE PLANAR ROTATOR 
WITH A GAUSSIAN APPROXIMATION 

By taking the second moment of Eq.(|7J), one obtains 
the following equation 



L/2 



-L/2 



dX I v I duj0(ujX - v)\v - LoX\F(u)(f>B{v)Auj 2 = 

(CI) 



where Aw 2 = uj* 2 — uj 2 . This equation means that for a 
stationary state the second moment of the distribution 
is time independent, or in other words, that the loss of 
the rotational energy of the planar rotator induced by 
inelastic collisions is compensated on average by collisions 
with bath particles with higher velocities. 
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By using Eq. JBJ , the difference between the square an- 
gular velocities at a collision is given by 



Aco 2 = -\(l + a) j 



^ + A 2 

m 

+ A 2 



\ m ) 



(C2) 



Introducing the dimcnsionless vectors 

m / / 



and 



2T V ' V 2T U 



(C3) 



G —[G X , Gy) 




(C4) 



Therefore the scalar product G.s gives 
G.s = v — lo\ 
Ea. ljCip can be expressed as 

rL/2 
'-L/2 



(C5) 



/■L/2 /■ 

/ dA / dscxp(-s 2 )6»(-G.s)|G.s| 

J-L/2 J 



2A(l + a)G.s /2T A 2 (l + a) 2 (G.s) 2 



^ + A 2 



V 771 



= 
(C6) 



Let us define a new coordinate system |21j where the y- 
axis is parralcl to G: units vectors are denoted (ei,e2) 
whereas the unit vectors of the original system (e x ,e v ). 
It follows that G = |G|ei and therefore one can write 
that 



Geie„ = IGlea; = G x 



(C7) 



Inserting Eq. (|C7|1 in Eq. I|C6(I allows for performing stan- 
dard gaussian integrals and finally one obtains Eq. 1311) 
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